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We report a study of disorder effects on epitaxial graphene in the vicinity of the Dirac point by 
magneto-transport. Hall effect measurements show that the carrier density increases quadratically 
with temperature, in good agreement with theoretical predictions which take into account intrinsic 
thermal excitation combined with electron-hole puddles induced by charged impurities. We deduce 
disorder strengths in the range 10.2 ~ 31.2 meV, depending on the sample treatment. We investigate 
the scattering mechanisms and estimate the impurity density to be 3.0 ~ 9.1 x 10^° cm“^ for 
our samples. An asymmetry in the electron/hole scattering is observed and is consistent with 
theoretical calculations for graphene on SiC substrates. We also show that the minimum conductivity 
increases with increasing disorder potential, in good agreement with quantum-mechanical numerical 
calculations. 

PACS numbers: 72.80.Vp, 71.23.-k, 72.10.-d 


Many of the exceptional electronic properties of 
graphene arise from its linear dispersion relation [1, 2]. 
However, when the Fermi energy approaches the Dirac 
point, its properties can be dominated by the effects 
of disorder, which can be both intrinsic (such as rip¬ 
ples, topological lattice defects) and extrinsic (including 
cracks/voids, adatoms, charged impurities, etc.), in gen¬ 
eral varying from sample to sample [3]. Of particular 
significance are the effects of disorder potentials on elec¬ 
trical transport properties [4] due to the lack of screening 
at very low carrier densities. Microscopically, the fluctu¬ 
ating electrostatic potential breaks up the intrinsically 
homogeneous charge distribution into electron-hole pud¬ 
dles [5-9]. This effect is recognised to mainly originate 
from unintendedly introduced charged impurities, whose 
type, spatial distribution and density also depend on the 
sample environment, device fabrication techniques, and 
particularly graphene synthesis and treatment processes. 

Recently, epitaxial graphene on SiC (SiC/G) has been 
reported to have very high quantum Hall breakdown 
current density [10] which potentially allows a quantum 
electrical resistance standard operating at even higher 
temperatures and lower magnetic fields [11]. Low and 
well controlled carrier density is required to achieve high 
breakdown current in these conditions, and understand¬ 
ing the disorder effects is therefore highly important. To 
date, there are very few experimental studies of disorder 
in epitaxial graphene grown on SiC due to the intrin¬ 
sically high level of doping from the substrate [12]. In 
this letter, using extremely low carrier density epitaxial 
graphene, we describe the role of disorder in governing 
the temperature dependent magneto-transport. 


Our SiC/G samples were epitaxially grown on the Si- 
terminated face of 4H-SiC. The devices used in this study 
all have an 8-leg Hall bar geometry of various sizes fabri¬ 
cated using standard electron beam lithography followed 
by O 2 plasma etching and large-area titanium-gold con¬ 
tacting. A non-volatile polymer gating technique was 
used to control the carrier density in epitaxial graphene 
by room-temperature UV illumination [13] or corona dis¬ 
charge [14]. Both DC and AC magneto-transport mea¬ 
surements were carried out using an Oxford Instruments 
21 T superconducting magnet with a variable temper¬ 
ature insert which allows temperature-dependent mea¬ 
surements from 1.4 K up to 300 K. 

Magneto-transport measurements were made on three 
SiC/G devices, which we denote GDI, CD2 and UVl. 
We used two different techniques to reduce the relatively 
high initial electron density and tune the Fermi level 
to the vicinity of the Dirac point: CDl and CD2 were 
treated with multiple negative ion projections produced 
by corona discharge onto the bilayer polymer gate, result¬ 
ing in extremely low final electron densities of 1.2 and 1.3 
X 10^^ cm“^, respectively; UVl was treated with UV illu¬ 
mination which eventually reduced the electron density 
to 8 X 10^^ cm“^. As we will show below, these val¬ 
ues should not be treated as the real electron densities, 
but merely are effective carrier densities, u-e//, calculated 
from the low-field Hall coefficients at 1.4 K assuming a 
homogeneous landscape with a single type of charge car¬ 
riers. Even so, these carrier densities should provide a 
good, rough estimation of the Fermi energy, which is be¬ 
tween 12.7 meV and 32.9 meV, based on the assumption 
of a linear dispersion where the density of states van- 
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FIG. 1. (Color online) The Hall resistance Rxy and the longitudinal resistance Rxx as a function of magnetic field at temperatures 
from 1.4 K to 300 K for sample CD2. The sample enters the quantum Hall regime already from about 0.6 T as observed from 
the quantised Rxy and the vanishing Rxx at low temperatures. 


ishes at the Dirac point [15]. In reality, due to the effects 
of disorder, a residual density of states and coexistence 
of electrons and holes [6] at Ep 0 are expected, thus, 
the determination of an extremely low Fermi energy from 
Hall effect measurements becomes non-trivial. Neverthe¬ 
less, it is evident that the actual Fermi energy should be 
lower than the above calculated values. The extremely 
low effective carrier densities and estimates of the Fermi 
level show that all measured devices can be treated as 
being effectively at the Dirac point. 

In Fig. 1 we present typical experimental results: the 
Hall resistance Rxy and the longitudinal resistance Rxx 
of sample CD2 as a function of magnetic field at temper¬ 
atures from 1.4 K to 300 K. In our study, all three devices 
show similar behaviour as shown in Fig. 1. Due to the 
extremely low carrier densities of the samples, quantum 
Hall plateau corresponding to the filling factor u = 2 
can be observed already from about 0.6 T at 1.4 K. The 
Hall resistance becomes significantly non-linear when ap¬ 
proaching the quantum Hall regime. Therefore, to ex¬ 
tract the zero-field carrier densities of our devices, only 
Hall coefficients between -0.1 T and +0.1 T are used. 

It has been theoretically studied and experimentally 
confirmed that, close to the Dirac point, as a consequence 
of disorder, the carrier density landscape is extremely in¬ 
homogeneous and electron-hole puddles are formed [4-9] . 
Classically, the low-field Hall coefficient in the presence 
of both electrons and holes is given by. 


When the Fermi energy is zero, i.e. at charge neutrality 
point (CNF), Ue = Uh > 0. Thus, Ug// = aue^ where 

^ • Notably, for electron-like behaviour [Rh < 

0), Of > 0; for hole-like behaviour {Rh > 0), a < 0. 

We now analyse the temperature dependence of the 
effective carrier density Ug// as shown in Fig. 2 for the 
three devices. A quadratic increase of rieff with increas¬ 
ing temperature can be clearly observed for all of the 
samples. Each sample also exhibits a distinct non-zero 
residual charge density at the low temperature limit even 
when Ef ^ 0, indicating that the potential landscape of 
our devices is highly inhomogeneous. These features are 
clearly different from the Arrhenius behaviour of conven¬ 
tional semiconductors and intrinsic thermal activation in 
graphene when no disorder effects are accounted for (i.e., 
there is no residual carrier density). Accurate fitting 
can be made based on the theory [4] assuming that the 
electronic potential energy of disordered graphene follows 
Gaussian statistics, which give the probability of finding 
the local potential within a range dV about E, 

P(V)dV= ^ p-^dV. (3) 

y2'KS^ 

where 5 is a parameter used to characterise the strength 
of the potential fluctuations. As a consequence, the 
temperature-dependent charge density at CNP for both 
electrons and holes are [4], 


^ 1 - nhPh 

~ FB e {uePe + nnyh)'^ ’ 


( 1 ) 


ne(T) = nh{T) 


gs9v {TThsT)'^ 

2Tr{hvF)^^4. 12 


( 4 ) 


where Ug (uh) and /ig {fih) are the electron (hole) den¬ 
sity and mobility, respectively. The carrier density di¬ 
rectly extracted from the low-field Hall effect is therefore 
effectively. 


where Qs = Qv = are the spin and valley degeneracies, 
and vp ~ 10^ m/s is the Fermi velocity. The temperature 
dependence of the effective carrier density is therefore, 

= <ang(T), (5) 


(fleRe ffi '^hRh) 
UeRl - UhrI 


where a is assumed to be constant over the tempera¬ 
ture range under consideration. The predicted temper- 
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TABLE L Energy Eluctuations of e-h Puddles in Graphene 


Synthesis (Treatment) 

Disorder Strength 

Probing Method 

Epitaxial on SiC (CDl) 

12.7 ± 0.6 meV 

Magneto-transport 

Epitaxial on SiC (CD2) 

10.2 ± 0.4 meV 

Magneto-transport 

Epitaxial on SiC (UVl) 

31.3 ± 2.0 meV 

Magneto-transport 

Epitaxial on SiC (AO) 

15 ± 1 meV 

Magneto-transport 

Epitaxial on SiC 

12 meV 

KPM [5] 

Exfoliated on Si02/Si 

50 meV 

SET [6] 

Exfoliated on Si02/Si 

~ 20 meV 

STM [7] 

Exfoliated on h-BN 

5.4 meV 

STM [8] 

CVD on Ir(lll) 

~ 30 meV 

STM/STS [9] 


EIG. 2. (Golor online) Temperature dependence of the ef¬ 
fective carrier densities Ueff deduced using Eq. (1) and (2) 
for sample GDI, GD2 and UVl. Quadratic increase with in¬ 
creasing temperature is observed, together with non-vanishing 
carrier densities ne//,o at T ^ 0 K. The experimental data 
is well fitted using Eq. (4) and (5) as shown in the figure 
(dash lines), where the disorder potential strength s and the 
mobility ratio [lej[ih are extracted from the fitting. 


ature dependence from Eq. (4) and (5) fits the experi¬ 
mental data very well (Fig. 2), giving potential fluctu¬ 
ation strengths s = 12.7, 10.2 and 31.3 meV, and pre¬ 
factors a which translate into mobility ratios of electrons 
to holes \iel— 2.04, 1.85 and 2.36, for the devices 
GDI, CD2 and UVl respectively. Table I shows compar¬ 
isons of the potential fluctuations due to electron-hole 
puddles, between the values deduced from our magneto¬ 
transport measurements and those found in the litera¬ 
ture [5-9] , where most of the characterizations are based 
on scanning tunnelling microscopy (STM). Table I also 
includes the disorder strength (15 ± 1 meV) from our 
analysis of the published data for SiC/G samples exposed 
to aqueous-ozone (AO) processing [16], which results in 
high mobility and extremely low p-type doping with an 
effective carrier density ne//,o = —4.0 x 10^^ cm“^ (neg¬ 
ative sign for hole-like behaviour) from Hall measure¬ 
ments. We find that the disorder strengths measured in 
our samples are consistent with those reported previously 
for SiC/G, and are smaller than those of CVD and exfoli¬ 
ated samples on Si 02 , while they are slightly larger than 
that of exfoliated graphene on h-BN, which is an atom¬ 
ically smooth, dangling bonds free and lattice-matched 
substrate to support high quality graphene [17]. These 
comparisons suggest that SiC/G generally has very good 
quality and relatively small amounts of disorder, even 
though the actual characteristics are expected to vary 
from sample to sample and may also be sensitive to the 
sample treatment as seen from Table I. At the same time, 
it is demonstrated that magneto-transport measurement 
is an additional effective method to investigate the dis¬ 
order effects and characteristics in graphene. 

To evaluate the scattering mechanisms in our SiC/G 
samples, we now turn to examine the temperature de¬ 


pendence of the longitudinal conductivity and the 
electron mobility as shown in Fig. 3. Carrier mobilities 
of individual species are calculated classically based on, 

Oxx = eiUeUe + nhUh) = + l^h), ( 6 ) 

a 

via Eq. (5) and the value ^ deduced from a. It is 
observed that (Jxx{T) remains slowly varying with weak 
non-monotonic fluctuations for a large range of temper¬ 
atures. Similar behaviour has been reported for mono- 
layer graphene samples when Ep ~0 [18, 19], and this is 
clearly different from thermally activated conductivity in 
conventional gapped semiconductors and from phonon- 
limited behaviour in graphene, which will result in a T~^ 
or dependence [20, 21] at low or high temperatures 
respectively due to intravalley acoustic phonon scatter¬ 
ing. It should be pointed out that this temperature de¬ 
pendence of conductivity in our extremely low carrier 
density samples could be a combination of various con¬ 
tributions. It is believed that this weakly varying conduc¬ 
tivity is mainly due to the temperature-dependent carrier 
density as described above and the /i(T) dependence as 
we will discuss below. At the lowest temperatures there 
are also temperature dependent weak localization correc¬ 
tions, which can be seen from Fig. lb around B = 0 T 
but have been excluded in Fig. 3a. Fig. 3b shows the 
electron mobility as a function of temperature, as well as 
the mobility limits as a result of various scattering mech¬ 
anisms, including impurity scattering, scattering by lon¬ 
gitudinal acoustic (LA) phonons in graphene, and also by 
remote interfacial phonons (RIP) at the SiC/graphene in¬ 
terface [22, 23]. In the case of charged impurities, carrier 
mobility is inversely proportional to the impurity density 
'^imp [15, 24], 

flimp ~ 5 V' / 

Tlimp 


where Co is a constant. For LA phonon scattering [20], 


_ ehpsVp% 
TmeD^ksT' 


(8) 


with ps = 7.6 X 10^ kg/m^ the two-dimensional mass 
density, Vs = 1.7 x 10^ m/s the sound velocity, and 
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Da = 18 eV the acoustic deformation potential. The 
RIP limited mobility is given by [22, 23], 


Mfl/p = — [T( 


Ci 


net 




( 9 ) 


where Ci and Ei are electron-phonon coupling constants 
and phonon energies of the phonon modes under consid¬ 
eration. To fit our data, we first considered three phonon 
modes: two out-of-plane acoustic phonon modes in epi¬ 
taxial graphene {Ei = 70 meV and E 2 = 16 meV) [ 22 , 25] 
and a surface phonon mode of 4H-SiC (E '3 = 117 meV) 
[22, 23, 26]. However, due to their relative large phonon 
energies, none of these can yield a reasonable fit, which 
can only be obtained (Fig. 3b solid lines) when an addi¬ 
tional low-energy phonon mode {E/^ ^ 2 meV) is intro¬ 
duced. This is consistent with the previously reported 
results [ 22 , 23, 27, 28], and this low-frequency remote 
phonon mode has been recognised to originate from the 
interaction between graphene and the buffer layer, that 
they are oscillating out-of-phase parallel to each other. 

It can be seen from Fig. 3b that impurity scattering 
plays the most dominant role at low temperatures (< 
100 K), while the high-temperature mobility is probably 
limited by RIP scattering, since LA phonons make only 
a small contribution to the overall mobility for tempera¬ 
tures below 400 K. Using Cq ~ 5 x 10^^ V“^s“^ [24], the 
densities of charged impurities for our SiC/G samples are 
estimated to be 3.0 ~ 9.1 x 10^^ cm“^, which are 1 ^ 
2 orders of magnitude lower than that in typical exfoli¬ 
ated [30] and CVD grown [31, 32] graphene on Si 02 , but 
are comparable to that of h-BN supported graphene [33] , 
consistent with its high charge carrier mobility. Even 
though we restrict the above analysis to phonon and im¬ 
purity scattering, other possible scattering mechanisms 
exist, such as scattering due to ripples [34, 35] and very 
large defects [36]. Quantitative analysis of these mecha¬ 
nisms on our devices is rather difficult since systematic 
examination of the sample morphology is required and, 
on the other hand, the theoretical pictures are rather 
complicated and still contentious [37]. 

So far we have been able to identify that charge carrier 
scattering at low temperatures in our SiC/G is mainly 
due to impurities, in the classical regime where quan¬ 
tum corrections are suppressed by magnetic fields. It is 
these impurities which provide the same origin for gen¬ 
erating the electron-hole puddles at Ep 0. Further¬ 
more, we believe that these impurities are most likely to 
be charged/Goulomb impurities rather than short-range 
impurities. The main evidence for this is the presence of 
unequal electron and hole mobilities, which is a conse¬ 
quence of the unbalanced transport cross sections due 
to charged scatterers in a system with 2D relativistic 
dispersion [24, 38]. This can be intuitively understood 
as because an attractive potential scatters a charge car¬ 
rier more effectively than a repulsive potential [38]. As 
presented in Fig. 2, we have obtained similar fie/l^h 
in the range of 1.85 ^ 2.36. According to the theory 
[38], assuming a single species of monovalent {\Z\ = 1) 
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FIG. 3. (Color online) (a) The longitudinal conductivity as a 
function of temperature, where weak non-monotonic depen¬ 
dences are shown, (b) The temperature dependence of the 
electron mobility of our samples. Individual contributions 
due to impurity scattering (green/pink/blue dash lines) for 
all three samples, LA phonon scattering (blue dash-dot line), 
RIP scattering (blue dash-dot-dot line) for UVl as an exam¬ 
ple are shown. The solid lines represent the overall /ie(T) 
dependence by fitting the experimental data, (c) cFmin as 
a function of disorder strength. An /3(s — A) 2 dependence 
(green dash-dot line) is observed from our experimental data 
(green triangles). <Jmin as a function of Kq from numerical 
calculations by Adam et al. (gray circles and solid line), as 
well as predictions using the Boltzmann (blue dash line) and 
the self-consistent Boltzmann (red solid line) theories are also 
shown [29]. 


impurities, the above mobility ratios can be translated 
into a dimensionless asymmetry factor c = 0.30 ^ 0.39, 
which is used to characterise the strength of this asym¬ 
metry effect (i.e. c = 0 for //g = and c ^ 1 for 
Me(/i) ^ M/i(e))- The nature of this asymmetry factor 
depends on the dielectric constant of the substrate: for 
Si 02 , c|e^= 3.9 ^ 0.46; for SiG substrates, the same as 
used in our devices, c|e^=io.o ~ 0.32, which is in very 
good agreement with our experimental results. Small 
variations around the predicted value are expected, since 
the actual electrostatic environment of each SiG/G sam¬ 
ple could also be affected by the polymer top-gate di- 
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electrics, meanwhile, the types and amounts of charged 
impurities present in our samples could be more complex. 

Finally, the effects of disorder potential fluctuations 
on the low-temperature non-vanishing minimum conduc¬ 
tivity (cTmin) at the Dirac point are investigated for 
graphene in the diffusive transport regime. This property 
has been extensively considered theoretically and the two 
main existing approaches lead to contradictory results 
[29]. The semiclassical Boltzmann transport theory pre¬ 
dicts a decreasing amin with increasing disorder strength. 
With a self-consistent modification to the Boltzmann the¬ 
ory, a subsequent increase of the minimum conductiv¬ 
ity for higher disorder strengths is predicted. On the 
other hand, the minimum conductivity treated quantum- 
mechanically [29, 39-41] is increased for the entire disor¬ 
der strength range for a non-interacting model using a 
Gaussian correlated disorder potential. Experimentally, 
very few studies can be found addressing this problem in 
the literature [24]. Shown in Fig. 3c is the minimum con¬ 
ductivity (at 5 = 0) as a function of the disorder strength 
s obtained from our measurements when quantum correc¬ 
tions have been taken into account, as well as theoretical 
predictions including the numerical calculation via the 
quantum mechanical approach by Adam et al [29], and 
results form the (self-consistent) Boltzmann theories, for 
L = 50^, where L is the sample length, ^ is the cor¬ 
relation length of the assumed random Gaussian poten¬ 
tial U{r) in the system and the dimensionless parameter 
Kq oc {U{r)U{r')) is the disorder strength used in the 
theories. Our experimental results show that the min¬ 
imum conductivity increases with increasing 5, roughly 
following a ^{s — A) 2 dependence locally in the (0.5 ^ 
2.5) X ^ range, highlighted by the green dash-dot line 
in the figure, where /3 and A are constants. This increase 
agrees qualitatively well with the theoretical predictions 
[29] from the quantum-mechanical approach, where we 
assume s oc \/Kq. However, our data do not agree with 


the results from the Boltzmann and the self-consistent 
Boltzmann theory as shown in the figure. In addition, we 
note that the minimum conductivity may have a complex 
dependence on the sample length and details of quantum 
interference effects [29, 42, 43], and also be a function of 
the charged impurity density riimp indicated from pre¬ 
vious experimental work by Ghen et al. [24], whose re¬ 
sults suggest that dmin drops with increasing riimp at low 
impurity densities and may saturate rapidly. To allow 
a more conclusive interpretation, however, more exper¬ 
imental data and systematic comparisons between well- 
controlled samples from different synthesis methods and 
a larger range of disorder potentials and impurity densi¬ 
ties would be needed. 

In summary, we have presented temperature de¬ 
pendent magneto-transport measurements on epitaxial 
graphene. We have demonstrated the disorder effects 
when the Fermi energy lies in the vicinity of the Dirac 
point, and have been able to identify the main origin 
of those effects to be charged impurities. The disorder 
strength and the impurity densities of our samples have 
been estimated from experimental results. We have also 
shown that the minimum conductivity increases with in¬ 
creasing disorder strength, in good agreement with nu¬ 
merical quantum-mechanical calculations. Overall, the 
application of this method can, therefore, provide an al¬ 
ternative and effective route for quantitatively studying 
the disorder characteristics in graphene and other two- 
dimensional materials. 

This work was partially supported by the U.K. EP- 
SRG and NMS, E.U. Graphene Flagship (Gontract No. 
GNEGT-IGT-604391), and EMRP GraphOhm. Work at 
Ghalmers University of Technology was supported by the 
Swedish Foundation for Strategic Research (SSF), Lin¬ 
naeus Gentre for Quantum Engineering, Knut and Alice 
Wallenberg Foundation and Ghalmers AoA Nano. J. H. 
acknowledges financial support from the Ghina Scholar¬ 
ship Gouncil. 


* r.nicholas@physics.ox.ac.uk 

[1] P. R. Wallace, Phys. Rev. 71, 622 (1947). 

[2] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 
M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and 
A. A. Firsov, Nature 438, 197 (2005). 

[3] A. H. Castro Neto, F. Guinea, N. M. R. 

Peres, K. S. Novoselov, and A. K. Geim, 
Rev. Mod. Phys. 81, 109 (2009). 

[4] Q. Li, E. H. Hwang, and S. Das Sarma, 

Phys. Rev. B 84, 115442 (2011). 

[5] A. E. Curtin, M. S. Fuhrer, J. L. Tedesco, R. L. 
Myers-Ward, C. R. Eddy Jr., and D. K. Gaskill, 
Appl. Phys. Lett. 98, 243111 (2011). 

[6] J. Martin, N. Akerman, G. Ulbricht, T. Lohmann, 
J. H. Smet, K. Von Klitzing, and A. Yacoby, 
Nat. Phys. 4, 144 (2008). 

[7] Y. Zhang, V. W. Brar, C. Girit, A. Zettl, and M. F. 
Crommie, Nat. Phys. 5, 722 (2009). 


[8] J. Xue, J. Sanchez-Yamagishi, D. Bulmash, P. Jacquod, 
A. Deshpande, K. Watanabe, T. Taniguchi, P. Jarillo- 
Herrero, and B. J. LeRoy, Nat. Mater. 10, 282 (2011). 

[9] S. C. Martin, S. Samaddar, B. Sacepe, A. Ki- 

mouche, J. Coraux, F. Fuchs, B. Grevin, 
H. Court ois, and C. B. Winkelmann, 

arXiv:1304.1183v3 [cond-mat.mes-hall] (2014). 

[10] J. A. Alexander-Webber, A. M. R. Baker, T. J. B. M. 
Janssen, A. Tzalenchuk, S. Lara-Avila, S. Kubatkin, 
R. Yakimova, B. A. Piot, D. K. Maude, and R. J. 
Nicholas, Phys. Rev. Lett. Ill, 096601 (2013). 

[11] A. Tzalenchuk, S. Lara-Avila, A. Kalaboukhov, S. Pao- 
lillo, M. Syvajarvi, R. Yakimova, O. Kazakova, 
T. J. B. M. Janssen, V. Fal’ko, and S. Kubatkin, 
Nat. Nanotechnol. 5, 186 (2010). 

[12] C. Riedl, C. Coletti, and U. Starke, 

J. Phys. D: Appl. Phys. 43, 374009 (2010). 

[13] S. Lara-Avila, K. Moth-Poulsen, R. Yakimova, 




6 


T. Bj0rnholm, V. Fal’ko, A. Tzalenchuk, and S. Ku- 
batkin, Adv. Mater. 23, 878 (2011). 

[14] A. Lartsev, T. Yager, T. Bergsten, A. Tzalenchuk, T. J. 
B. M. Janssen, R. Yakimova, S. Lara-Avila, and S. Ku- 
batkin, Appl. Phys. Lett. 105, 063106 (2014). 

[15] E. H. Hwang, S. Adam, and S. D. Sarma, 
Phys. Rev. Lett. 98, 186806 (2007). 

[16] T. Yager, M. J. Webb, H. Grennberg, R. Yaki¬ 
mova, S. Lara-Avila, and S. Kubatkin, 

Appl. Phys. Lett. 106, 063503 (2015). 

[17] C. R. Dean, A. F. Young, 1. Meric, C. Lee, L. Wang, 
S. Sorgenfrei, K. Watanabe, T. Taniguchi, P. Kim, K. L. 
Shepard, and J. Hone, Nat. Nanotechnol. 5, 722 (2010). 

[18] K. 1. Bolotin, K. J. Sikes, J. Hone, H. L. Stormer, and 
P. Kim, Phys. Rev. Lett. 101, 096802 (2008). 

[19] S. V. Morozov, K. S. Novoselov, M. 1. Katsnelson, 
F. Schedin, D. C. Elias, J. A. Jaszczak, and A. K. Geim, 
Phys. Rev. Lett. 100, 016602 (2008). 

[20] E. H. Hwang and S. Das Sarma, 
Phys. Rev. B 77, 115449 (2008). 

[21] T. Stauber, N. M. R. Peres, and F. Guinea, 

Phys. Rev. B 76, 205423 (2007). 

[22] A. J. M. Giesbers, P. Prochazka, and G. F. J. Flipse, 
Phys. Rev. B 87, 195405 (2013). 

[23] S. Tanabe, Y. Sekine, H. Kageshima, M. Nagase, and 
H. Hibino, Phys. Rev. B 84, 115458 (2011). 

[24] J.-H. Ghen, G. Jang, S. Adam, M. S. Fuhrer, E. D. 
Williams, and M. Ishigami, Nat. Phys. 4, 377 (2008). 

[25] J. Cervenka, K. van de Ruit, and G. F. J. Flipse, 
Phys. Status Solidi B 247, 2992 (2010). 

[26] H. Nienhaus, T. U. Kampen, and W. Monch, 

Surf. Sci. Lett. 324, L328 (1995). 

[27] G. Yu, J. Li, Q. B. Liu, S. B. Dun, Z. Z. 

He, X. W. Zhang, S. J. Gai, and Z. H. Feng, 


Appl. Phys. Lett. 102 , 013107 (2013). 

[28] N. Ray, S. Shallcross, S. Hensel, and O. Pankratov, 
Phys. Rev. B 86, 125426 (2012). 

[29] S. Adam, P. W. Brouwer, and S. Das Sarma, 
Phys. Rev. B 79, 201404 (2009). 

[30] Y.-W. Tan, Y. Zhang, K. Bolotin, Y. Zhao, S. Adam, 
E. H. Hwang, S. Das Sarma, H. L. Stormer, and P. Kim, 
Phys. Rev. Lett. 99, 246803 (2007). 

[31] J. Lee, L. Tao, Y. Hao, R. S. Ruoff, and D. Akinwande, 
Appl. Phys. Lett. 100, 152104 (2012). 

[32] J. Ghan, A. Venugopal, A. Pirkle, S. McDonnell, D. Hi- 
nojos, G. W. Magnuson, R. S. Ruoff, L. Golombo, R. M. 
Wallace, and E. M. Vogel, AGS Nano 6, 3224 (2012). 

[33] K. M. Burson, W. G. Gullen, S. Adam, G. R. Dean, 
K. Watanabe, T. Taniguchi, P. Kim, and M. S. Fuhrer, 
Nano Lett. 13, 3576 (2013). 

[34] F. Guinea, J. Low Temp. Phys. 153, 359 (2008). 

[35] M. 1. Katsnelson and A. K. Geim, 
Phil. Trans. R. Soc. A 366, 195 (2008). 

[36] J.-H. Ghen, W. G. Gullen, G. Jang, M. S. Fuhrer, and 
E. D. Williams, Phys. Rev. Lett. 102 , 236805 (2009). 

[37] S. Das Sarma, S. Adam, E. H. Hwang, and E. Rossi, 
Rev. Mod. Phys. 83, 407 (2011). 

[38] D. S. Novikov, Appl. Phys. Lett. 91, 102102 (2007). 

[39] H. Suzuura and T. Ando, 

Phys. Rev. Lett. 89, 266603 (2002). 

[40] A. Schuessler, P. M. Ostrovsky, 1. V. Gornyi, and A. D. 
Mirlin, Phys. Rev. B 79, 075405 (2009). 

[41] J. Tworzydlo, G. W. Groth, and G. W. J. Beenakker, 
Phys. Rev. B 78, 235438 (2008). 

[42] F. Ortmann, A. Gresti, G. Montambaux, and S. Roche, 
Europhys. Lett. 94, 47006 (2011). 

[43] N. Leconte, A. Lherbier, F. Varchon, P. Or- 

dejon, S. Roche, and J.-G. Gharlier, 

Phys. Rev. B 84, 235420 (2011). 


